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Abstract 

A simple technique for expanding the free energy of general six- vertex 
models about free-fermion points is introduced. This technique is used 
to verify a Coulomb gas prediction about the behavior of the leading 
singularity in the free energy of the staggered F-model at zero staggered 
field. 



1 Definition of the staggered F-model 

The staggered F-model is a special case of the six- vertex model. The six- vertex 
model can be defined as follows: place arrows on the edges of a square lattice 
so that there are two arrows pointing into each vertex. Six types of vertices 
can arise (hence the name of the model). These vertices are shown in fig. ^ 
By giving each vertex-type a (position-dependent) energy the model is defined. 
These models were first introduced to study (anti-)ferroelectric systems. Later 
it was shown that six-vertex models can be mapped to solid-on-solid models 
13] ■ Only a few of these models can be solved exactly. These include the free- 
fermion models and models that can be solved using the Bethe Ansatz 
El Ell EH 12] • To define the staggered F-model, we divide the lattice into two 
sublattices A and B, such that the nearest neighbor of an A vertex is a B vertex. 
The vertex energies are chosen as indicated in fig. When the staggered field 
(s) vanishes the model reduces to the F-model, which has been solved by Lieb 
[§]• At zero staggered field the model is critical. In this case the groundstate is 
twofold degenerate consisting of vertices of type 5 on sublattice A and vertices 
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Figure 1: The six vertices and their energies. The upper and lower signs correspond 
to sublattice A resp. sublattice B. 
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of type 6 on sublattice B, or vice versa. For f3e > a nonzero staggered field 
lifts this degeneracy, and forces the model into an ordered state |I5| . 



2 Coulomb gas results 

By assuming that the F-model renormalizes to the Gaussian model, it is possible 
to find the behavior of the staggered F-model in infinitesimal staggered fields 
It is found that the leading singularity in the free energy is: 



F s (f3e,/3s) « {(3s) 2 -^m (i) 

where 

j 09e) = \ arccos (l - ^ exp (2/?e)J (2) 

At the point fie = \ In (2 — \/2) ~ —0.2674 the exponent becomes infinite. 
Below this point a finite staggered field is necessary to force the model to an 
ordered state. In this case the transition to the ordered state happens via a 
Kosterlitz-Thouless (KT) transition. The existence of a line of KT transitions 
intersecting the point (fie — \ In (2 — \/2) ,0s = O) has been verified by com- 
bining the results of transfer matrix studies with scaling arguments 

In this paper we will verify Q by expanding around Baxter's exact solution 
on the line f3e = \ In 2 



3 Baxter's solution of the staggered F-model 

Baxter has solved the staggered F-model at fie = | ln(2) 1 . Later it was found 
that this solution could be generalized to other models if a certain condition 
involving the vertex weights is met. This condition is called the free-fermion 
condition because for eight-vertex models satisfying this condition the problem 
leads to a problem of noninteracting fermions in the S-matrix formulation. Let 
u>i be the vertex weight for a vertex of type i (see fig. ^1 , then the free-fermion 
condition for six- vertex models is: 

w\W2 + W3W4 — W5W6 = (3) 

The weights Wi may be chosen inhomogeneous. We now proceed by presenting 
Baxter's solution of the staggered F-model. 

Divide the lattice into two sublattices A and B. Choose the vertex energies as 
indicated in fig. ^ Consider the ground state in which all A- vertices are vertices 
of type 6, and all B-vertices are of type 5. Any state can now be represented 
by drawing lines on the lattice where the arrows point oppositely to the ground 
state configuration. In terms of these lines the six vertices are represented by 
vertices with either no lines, two lines at right angles, or four lines. The energies 
of these vertices are respectively — s, e and s. The next step is to replace the 
original lattice by a decorated lattice by replacing each original vertex by a 
"city" of four internally connected points (see fig-EJ. 

The lines on the original lattice are regarded as dimers on the external edges 
of the decorated lattice. For any configuration on the original lattice, it is 
possible to place dimers on the internal edges of the decorated lattice, so that 
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Figure 2: The "cities" on the decorated lattice. A and B refer to the two sublattices. 
The meaning of the orientations on the edges is explained in the text. 



the lattice becomes completely covered. Now associate to each dimer a weight as 
indicated in fig. El Demanding that the closed-packed dimer problem formulated 
on the decorated lattice is equivalent to our original problem yields: 

C = expf-i/fr) (4) 

u = iV2ex P Q/3^ (5) 

pe = -In (2) (6) 

To solve the close-packed dimer problem, we use the Pfaffian method [HI 
I12| . This method expresses the partition function Z for a closed packed dimer 
model on an N by M planar lattice: 

Z 2 = det R (7) 

Here R is an N X M by N X M anti-symmetric matrix, defined as follows. 
Enumerate all the N x M vertices on the decorated lattice. If vertex i is not 
connected to vertex j via an edge, Rij = 0, else Rij = ± fugacity of dimer at 
edge connecting i to j. The way the signs have to be chosen is explained 0. 
These signs define an orientation of the edges. Positive Rij is indicated by an 
arrow pointing from i to j. 

To set up a perturbation theory about (3e = |ln(2), we also need the in- 
verse of R. Both the determinant and the inverse of R are easily calculated 
by performing a similarity transformation, see £Q for details. The determinant 
yields the following expression for the reduced free energy per vertex (i.e. the 
free energy times — /3), denoted as ^Baxter, for an infinite by infinite lattice: 



^Baxter = limA^Af^oo 2 N M l n det R 

fill r 27T , rr, 1 fnr, \ , r, / /, \ In M 7/1 m \° j 



+ 2NM 

W Jo* Jo* ln I 2 cosh ( 2 ' 3s ) + 2 cos (°i) cos d9id9 2 
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4 Perturbation theory 



We now proceed with the derivation of a perturbation theory about the free- 
fermion line of a six-vertex model. The Hamiltonian of a general six-vertex 
model can be defined as follows. One assigns an energy e (p, i) to a vertex in 
state p (see fig.nj and position i. The configuration of the lattice can be specified 
by a function c which maps a position of a vertex to a number, 1 • • • 6, which 
is to be interpreted as the state of the vertex at that position. The reduced 
Hamiltonian (H) is defined to be the functional that assigns to each state c its 
energy times —j3. We can thus write 

£r(c) = -/?5>( c (*),i) (9) 

i 

For H a Hamiltonian of a general six-vertex model and Hq a Hamiltonian 
of a free-fermion model, a perturbation V can be defined so that we have 

H = H Q + V (10) 

The partition function Z can be written as: 

Z = ex P (Ho (c) + V (c)) = Z (exp (V)) (11) 

c 

Here Zq is the partition function of the free-fermion model. The reduced free 
energy can be expressed as: 

F = Fq + In (exp (V)) = F Q + (V) + ^ (pT - (V)) 2 ) + ... (12) 

Here Fq is the reduced free energy of the free-fermion model. Now write V = 
Vi with Vi (c (i)) a perturbation of the vertex energy times — (5 at position i. 
(112(1 can be rewritten as: 

F = F Q + ]T (Vi) + \ [(ViVj) (Vi) (V,)} + ... (13) 

i " ij 

To compute a free-fermion average (Vi ± Vi 2 ...Vi n ), we can proceed as follows: 
Introduce a constraint in the free-fermion model by requiring the vertices at the 
positions i\ . . . i n to be in the states x\ . . . x n . The partition function of this 
model is denoted by Zi t ..^ n (x\ . . . x n ). We can then write 

<^ 2 ...VO= E ^-^■■■^l V ^-- V ^ (14) 

X\.. .x n 

It now remains to calculate Z^...^ (xi . . .x n ). It is convenient to reformulate 
this problem as follows: Denote the state of an arrow located at the edge j by 
Sj . Put Sj = 1 if the arrow points oppositely to the ground state configuration 
and Sj = otherwise. Define a constrained free-fermion model by requiring the 
arrow at the edge j r to be in state Sj r for 1 < r < m. We then want to evaluate 
the partition function of this model, which we denote as Z cons (sj 1 . . . Sj m ). The 
idea is to perturb the weights of the dimers on the edges j r infinitesimally. We 
redefine the weight of the dimer on the edge j r by multiplying it by (1 + e r ). 
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The partition function of the redefined free-fermion model (Z (ei . . . e m )) can be 
written in terms of the constrained partition functions as: 



Z(e 1 ... e m ) = Z c ° ns (s h ■ ■ ■ s jm ) JJ (1 + s jh e k ) 

{«} k=i ( 15 ) 

= Z + £ ^ cons = 1) e k + Z C ° nS (*A = 1. *ii = 1) ^ + . . . 



fc<; 



Z (ei . . . e,„) can be calculated using JJJ), by making the necessary changes to R. 
We can write: 



R — Rq + e kR( 



k) 



(16) 



fe=i 



Here Rq is the original unperturbed matrix, R^ is defined as follows: 

R(k)Sj — Ro,ij 

if i and j are connected by jk 

R(k),i 3 = 

if i and j are not connected by jk ■ 

Note that the R(k) have only two nonzero matrix elements. Inserting in 
(jTJ) and expanding gives: 



Z (ei . . . e m ) = Vdeti? = ^det i? exp I - Trln 



1 + X/ fci? o ^(fe) 



= \/ dct Rq 



1 + lzZ e * Tl ( R o lR W 



- Tr (R^Rw) Tr (iZ^JZp)) - Tr (R^R (k) R^R (l) ) 



(17) 

Using l|17|l and (|15|) we can directly read off the constrained partition functions 
if all the constrained arrows point oppositely to the ground state configuration. 
To calculate a general constrained partition function one can apply the principle 
of inclusion and exclusion. E.g. consider the evaluation of Z (si, s 2 , S3, S4, S5), 
with si = s 2 = 1 and S3 = S4 = S5 = 0. Put £3 = <4 = £5 = 1. According to the 
principle of incusion and exclusion, we can write: 

Z(si,s 2 ,s 3 ,S4, s 5 ) = Z(si,s 2 ) 

- [Z (s 1 ,S 2 ,t 3 ) +Z(8l,8 2 ,U) + Z (S!, s 2 ,t 5 )] 

+ Z {s 1 ,s 2 ,t 3 ,ti} + Z (si,a 2 !<3,*5) + Z (si,S2,t 4 ,t 5 ) 

- Z(si,s 2 ,t 3 ,U,t 5 ) (18) 



5 First order computation for the staggered F- 
model 

For the staggered F-model the expansion can be simplified. The vertex in the 
ground state at a particular point will be referred to as an a-vertex. A b-vertex 
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is obtained by reversing the arrows of an a-vertex. An a-vertex (b-vertex) is 
thus of type 5 or 6 and has an energy of — s (s). The constrained partition 
function corresponding to the model with one vertex constrained to be an a- 
vertex (b-vertex) is denoted as Z a (Z^). Note that under the transformation 
s — > — s the role of vertices a and b are interchanged. We thus have 



Z a (/3s) = Z h (-(3s) 



(19) 



If we put 0e = | In (2) + U we have, according to lll.ifl and 1)14(1. to first order 
in U: 

F = F - Z °~ Z *~ Zh U + (U 2 ) (20) 
Zo 

Here F is the reduced free energy per vertex of the staggered F-model, and 
Fq = ^Baxter in (JHJ ■ To calculate Z^ we only have to constrain two opposing 
arrows of one vertex to point oppositely to an a-vertex. Using the formalism of 
the previous section, we have obtained: 



Zo 
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27r r d9 lC ie 2 exp { ~ 2f3s) + cos {0l) cos {02) 

o Jo 



64tt 4 [7 J cosh (20s) + cos (0i) cos {6 2 . 

Using this, the first order expansion of the free energy can be written as: 
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F = F 



dFp 
80s 



1 



U 



(21) 



(22) 



6 Singular behavior in the vicinity of the free- 
fermion line 

We will now verify the Coulomb gas result (see section EJ: 



F s ~(0s) 2 ~^T (23) 

where 

j (0e) = I arccos f 1 - i exp (20e) J (24) 



If we put 

0e = X - In (2) + U (25) 
Expanding in powers of U yields: 



F S =A(U) (0s) 2 



-(U + O (U 2 )) In (0s) + ^(U 2 + (U 3 )) In 2 \0s\ + 

TT TT Z 



(26) 

where the amplitude A (U) is a meromorphic function. If we compare this with 
the non-analytical behavior at U — (see ()4()|l in the appendix), we find: 

A(U) = ±j + 0(l) (27) 

It then follows that the amplitude of the term (0s) 2 In 2 \0s\ is ^ (U + O (U 2 )). 
It is now a simple matter to verify this using (|22|l and l|40|) . From <|40ll and l|22|) 
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it follows that the order U contribution to the singular part of the reduced free 
energy, F% (/3s), can be written as 

F hs (/3s) = [Bi (/3s) In |/3s| + B 2 (/3s) In 2 |/3s|] U (28) 

with B\ and B 2 regular functions of /3s. Inserting (|40|l in (|22|l gives 



B 2 (/3s) 



2 7Q 
(/3s) 2 ~-(/3s) 4 + -(/3s) 6 + ... 



(29) 



3 v ' 90 
We have thus verified (|24|l to first order in U. 

7 Conclusions and outlook 

We have presented a simple technique for expanding the free energy of six- vertex 
models about free-fermion points. Applying this technique to the staggered F- 
model has enabled us to verify a Coulomb gas prediction about the singular 
part of the free energy of this model. It would be interesting to perform such 
computations to higher order in the free-fermion expansion. It is possible that 
such an undertaking might lead to proofs of certain Coulomb gas results. 

A Singular part of the free energy 

In this appendix we calculate the singular part of the free energy of the staggered 
F-model at /3s = on the free-fermion line. Expanding the logarithm in (jHJ 
yields 



iw vis) = dhdth s~ i cos ; (9i)cos " (92) 

(2n)\ 



2 2^n=l 2n cosh 2 ™ (2/?s) 



4"n! 2 



' = ^ CXP (-) ( E (3D 



Using the asymptotic expansion 

B 2 k 

2k 

\k=l 

where the B r are the Bernoulli numbers, we find 



oo 



1 

Baxter ' 



4tt ^ n 2 cosli {20s) 



1 1 1 



We can find the non-analytical part of the function Y^=i n p cosh*" (20s) as 



An 32n 2 128n 3 

l 

nP cosh 2 "(2/3s) 



(32) 



lows: Put t = In (cosh 2 (2/3s)) . We then have to find the non-analytical part of 
the function U p (t) with 

f,W-f=t!Si (33) 



n=l 



at t = for p > 2. From J33J| it follows that 



A = ~ U > (34) 
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We denote the non- analytical part of U p by U p . It then follows from 134( 1 that 

= -U P (35) 



dUp+i 



dt 

For p = 1 the sum in ((33(1 is easily evaluated: 

Ut{t) = - In (1 - exp (-*)) 
And we see that U\ (t) is given by 

Ux{i) = -In(t) 
From (j2Z| and JSlJl it then follows that 

^W = (- 1 ) P 7^4rr ln W 



(p-i)l 



Inserting this in ((32(1 gives 



— t H 1 

4tt V 8 192 3072 



m(f) 



(36) 
(37) 

(38) 
(39) 



Where F s (/3s) is the singular part of the free energy and t = 2 In (cosh (2/3s)). 
Expanding 1(39(1 in powers of f3s gives 



ft (0s) 



180 



5040 



In \/3s\ (40) 
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